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The great potential of Dirac electrons for plas- 
monics and photonics has been readily recog¬ 
nized after their discovery in graphene, followed 
by applications to smart optical devices. Dirac 
carriers are also found in topological insulators 
(TI) quantum systems having an insulating gap 
in the bulk and intrinsic Dirac metallic states 
at the surface . Here, we investigate the plas- 
monic response of ring structures patterned in 
Bi2Se3 TI films, which we investigate through ter¬ 
ahertz (THz) spectroscopy. The rings are ob¬ 
served to exhibit a bonding and an antibond¬ 
ing plasmon modes, which we tune in frequency 
by varying their diameter. We develop an an¬ 
alytical theory based on the THz conductiv¬ 
ity of unpatterned films, which accurately de¬ 
scribes the strong plasmon-phonon hybridization 
and Fano interference experimentally observed as 
the bonding plasmon is swiped across the promi- 
neng 2 THz phonon exhibited by this material. 
This work opens the road for the investigation of 
plasmons in topological insulators and for their 
application in tunable THz devices. 

Plasmons -the collective oscillations of charge carriers 
in conducting materials- hold great potential for com¬ 
bining electronics and photonics at the nanoscale. These 
excitations can propagate along extended surfaces or they 
can be localized at so-called hotspots. The engineer¬ 
ing of plasmon frequencies and spatial profiles has been 
mastered in metallic micro- and nanostructures such as 
nanospheres, nanorods [T], dimers mm, particle arrays 
0113 , and a plethora of more exotic morphologies [OHn]. 
Plasmon hybridization [ME], Fano resonances [HIM, 
and electromagnetically induced transparency |T3| are 
among the feats that have been realized and broadly used 
to understand and design plasmonic devices. The range 
of applications of plasmon excitations is vast and in¬ 
cludes optical sensing [ME], quantum electrodynamics 


[MUD], nonlinar optics EHE], photovoltaic technolo¬ 
gies [13], and medical diagnosis and treatment [331135] . 

Extensive experimental efforts are currently underway 
to find materials with improved plasmonic performance, 
in particular in the mid-infrared and terahertz parts of 
the electromagnetic spectrum. Examples of such ma¬ 
terials are low- m and high-Tc superconductors m, 
conductive oxides m, and graphene [351131] . The lat¬ 
ter exhibits a peculiar electronic structure, which en¬ 
ables unprecedented levels of electro-optical tunability 
via chemical or electrostatic doping [351137] : electrons 
in graphene behave as massless Dirac fermions charac¬ 
terized by a linear dispersion relation and a vanishing 
density of states at the Fermi energy, so that a few addi¬ 
tional charge carriers produce a large shift in the chemical 
potential. Dirac charge carriers are also found in three- 
dimensional topological insulator materials -i.e., quan¬ 
tum systems characterized by an insulating electronic gap 
in the bulk, which opening is due to strong spin-orbit in¬ 
teraction, and gapless surface states at the edge and in¬ 
terface. Surface states in TIs are metallic and chiral, as¬ 
sociated with massless Dirac quasi-particles [SHHID], 
protected from back-scattering by time-reversal symme¬ 
try. This means that they cannot be destroyed or gapped 
by scattering processes, other than magnetic impurities. 
Furthermore, unlike graphene, TI surface states sponta¬ 
neously provide a 2D Dirac-fermion system, segregated 
from the bulk material without the need of physically im¬ 
plementing an atomic monolayer. Since their discovery, 
TFs have attracted a growing interest due to their po¬ 
tential application in quantum computing [MIID], THz 
detectors [13], and spintronic devices [11], and some of 
these applications could benefit from propagating and 
localized plasmons sustainted by these materials. 

Recently, propagating plasmons have been observed 
in the THz domain in micro-ribbon arrays patterned 
on Bi 2 Se 3 TI films [15]. The plasmon frequency vp 
was found to follow an optical dispersion relation vp oc 
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1/vW, where W is the ribbon width, characteristic of 
2D metal systems and in excellent agreement with that 
predicted for Dirac plasmons. @51117! However, despite 
their potential for novel applications, neither localized 
plasmons nor their hybridization and coupling have yet 
been explored in TI micro and nano structures. In this 
Letter, we reveal the behavior of localized plasmons in 
topological insulator rings (TIRs), with a particular fo¬ 
cus on their dependence on size and temperature, as well 
as on their hybridization with a Bi2Se3 optical phonon 
located at 2 THz. We demonstrate geometrical tuning of 
the TIR plasmons and the resulting hybridization that 
exhibit a strong temperature dependence as observed in 
transmission spectra. This study indicates that TIRs 
hold a great potential for applications in thermally con¬ 
trolled THz devices. 


RESULTS 

Five films of Bi2Se3 (thickness h = 60 quintuple lay¬ 
ers ~ 60 nm) were grown by molecular beam epitaxy on 
0.5mm thick sapphire (AI2O3) substrates [HISS]- The 
films were characterized through resistivity and Hall mea¬ 
surements @H1 @5] (see Methods for details) indicating 
the presence of a 2D gas of Dirac electrons with a surface 
density nn = 3±1 x 10^^ cm“^. One of the films was kept 
as grown to serve as a reference, while the other four were 
patterned by electron-beam lithography and reactive-ion 
etching to form periodic hexagonal arrays of rings with 
different sizes. The outer (D) and inner (d) ring diam¬ 
eters were chosen to have a similar disk-to-hole aspect 
ratio (D/d ~ 3), while the lattice period a was taken to 
be approximately 4 times the inner diameter (a/d = 4) 
in all samples. Fig. shows a scanning electron mi¬ 
croscope (SEM) image of the D = 15 fim patterned film, 
whereas Fig. [^3 presents optical microscope images of 
the different arrays. 

The transmittance T(v) of the five films was measured 
in the terahertz domain using a Fourier-transform in¬ 
terferometer at temperatures held in the range of 10 to 
300 K in the AILES beamline at the synchrotron Soleil 
jsniisi]- Each transmittance has been normalized to that 
of the sapphire substrate collected at the same temper¬ 
ature. The corresponding extinction coefficient (Extinc¬ 
tion = 1-Transmittance) is plotted in Fig. for the 
as-grown sample and in Fig. for the patterned films 
at 10K (blue curves), 100K (green curves), and 300K 
(red curves). 

The extinction spectra in Fig. Bi exhibit several peaks 
shifting to higher frequency (blueshifting) as the ring 
size is decreased. These features are accompanied by 
a prominent dip-peak structure at a nearly constant fre¬ 
quency ~ 2 THz, independent of ring size. Additionally, 
a strong dependence on temperature is observed, includ¬ 
ing a reduction in the contrast of the observed dip-peak 


structure, as well as a broad increase in the extinction at 
higher frequencies. 

Quantitative insight into the origin of the observed 
spectra is obtained from theoretical simulations per¬ 
formed starting from the dielectric properties of the un¬ 
patterned film (see Methods for details). Since the thick¬ 
ness of the Bi2Se3 film is small compared with both the 
lateral size of the rings and the terahertz wavelengths, we 
can model them through a 2D surface conductivity (see 
Methods). Before analyzing the effects of nanostructur¬ 
ing, it is worth looking into the properties of this conduc¬ 
tivity, as extracted from extinction measurements for the 
unpatterned film of the same thickness (Fig. &)• The 
extinction of this film exhibits a phonon peak (a mode) 
at 1.85 THz, already observed in Bi2Se3 single crystals 
|52] . This phonon is superimposed on an electronic back¬ 
ground that decreases monotonically with frequency and 
that is associated with Dirac surface states @3]. The 
resulting conductivity is then modeled through two con¬ 
tributions: a Lorentzian to represent the a phonon and 
a Drude term accounting for the response of the surface 
Dirac electron gas (see Eq. B in Methods). The for¬ 
mer is expected to depend linearly on the film thickness, 
as the optical phonons are spatially localized and con¬ 
tribute all over the volume of the material. In contrast, 
the Drude term is independent of thickness but depends 
on the combined surface electron densities on both sides 
of the film. This metallic contribution is also strongly 
dependent on temperature (see Fig. presumably as 
a result of thermal smearing of the surface electron gas 
(notice that KbT/H = 2.1 THz at T = 100 K). 

By fitting the spectra of the homogeneous film (see 
Tabin Methods), we can parametrize the TI film con¬ 
ductivity and use it to perform classical electromagnetic 
simulations as discussed in the Methods section. Calcu¬ 
lations are carried out for all three temperatures under 
consideration and the resulting theoretical spectra very 
well reproduce both the experimental extinction of the 
unpatterned film (Fig. and those of the patterned 
ones (Fig. &)• In particular for the ring-arrays (Fig. 
analytical calculations reproduce the extinction shape, 
the main features of the temperature dependence noted 
above, as well as the general trend of the evolution of 
spectral maxima and minima with ring size. 

If we artificially switch off the a phonon contribu¬ 
tion from the film conductivity (i.e., by removing the 
Lorentzian term in Eq. B of Methods), the extinction 
spectra are simplified and consist of two dominant max¬ 
ima that evolve towards higher frequencies for smaller 
rings (Fig. [^, dashed curves). Because only the Drude 
term contributes to the conductivity in this case, the 
observed features can be undoubtedly assigned to plas- 
monic excitations. Closer examination of the origin of 
these features (see Methods) reveals that they are the 
bonding (at low energy) and antibonding (at higher en¬ 
ergy) ring plasmons. Moving back to the full simulation 
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FIG. 1: Geometrical tuning of the plasmonic response of topological insulator rings (TIRs). (a) Detailed SEM 
image of a Bi 2 Se 3 TIR having an outer (inner) diameter D= 15 (d = 5) ^m. (b) Extinction (= 1— Transmittance) spectra 
of unpatterned Bi 2 Se 3 film (60 nm thickness) at different temperatures. Dashed lines represent calculated extinctions for the 
unpatterned film (see text and Methods), (c-e) Optical images (c), measured extinction (d), and calculated extinction (e) of 
Bi 2 Se 3 ring arrays for different ring sizes, as specified in (c) (see labels with the values of D and d, respectively). The array 
lattice period is 4d and the D/d aspect ratio is close to 3 in all samples, whereas the ring size decreases from the top to the 
bottom panel. The transmittance is normalized to that of the sapphire substrate measured at the same temperature. The 
calculations of (e) are performed with (solid curves) and without (broken curves) the inclusion of the TI a phonon. 


including phonons (Fig. [^, solid curves), we still iden¬ 
tify the presence of both types of plasmons, although the 
broad bonding mode is now strongly perturbed by the rel¬ 
atively narrow optical phonon. The resulting plasmon- 
phonon interaction gives rise to characteristic dip-peak 
line shapes that are reminescent of the Fano interference 
observed in other plasmonic systems [Hilliisjisi]. 


For the smallest ring under consideration (Z? = 3 /rm), 
a slight disagreement is observed at high frequency be¬ 
tween experimental data and theory (cf. Fig. [^,e, 
bottom panels). In particular, the experimental anti¬ 
bonding absorption peak (above the a phonon), appears 
to be broader than the calculated one. This could be as¬ 
cribed to unaccounted decay channels other than ohmic 
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(single-particle) scattering, which is the only one consid¬ 
ered in the calculation and common to both plasmons. 
In this respect, fabrication defects cannot be ruled out, 
as they are expected to have stronger influence on the 
high-frequency (antibonding) plasmon mode. 


A broad picture of plasmon-phonon hybridization and 
its dependence on ring size is presented in Fig. where 
the background color plot shows the extinction coeffi¬ 
cient obtained from our analytical model (see Methods), 
including the TI a phonon, as a function of inverse outer 
ring diameter l/D and radiation frequency v. This is 
compared with the peak frequency of bonding (lower red- 
dashed curve) and antibonding (upper red-dashed curve) 
plasmons obtained when the a phonon is dismissed; the 
bonding mode traverses the a phonon frequency (black- 
dotted horizontal line) and these two modes hybridize, 
giving rise to extinction maxima above and below the 
phonon line (and a minimum around the bare phonon 
frequency). In contrast, the antibonding plasmon is far 
up in frequency compared to the phonon, so that it only 
experiences a mild redshift due to interaction with the 
latter. The experimentally observed maxima and min¬ 
ima (Fig. solid and open symbols, respectively) are in 
excellent agreement with the analytical calculation (color 
plot), including a similar hybridization of the bonding 
plasmon and the phonon, and a lack of hybridization of 
the antibonding plasmon. This bonding plasmon-phonon 
hybridization is a classical example of avoided crossing, 
which can be actually described through a simple inter¬ 
action model using an effective Hamiltonian of matrix 
elements (pl|i?|pl) = 27r/ii/pi, (ph|iJ|ph) = 2nhvpi^, and 
(pl|iF|ph) = (ph|7J|pl) = 2nhi'int, where |pl) and |ph) 
represent the plasmon and phonon modes, respectively, 
and t'int describes their interaction. Taking the plasmon 
and phonon frequencies Vp\ and from the lower red- 
dashed curve and the black-dotted line of Fig. re¬ 
spectively, we find the following hybridized frequencies 
v± = (1/2) Vpi + Vph±yJ (r-pi - izph)^ + . These 

are shown as solid (red) curves in the same figure giv¬ 
ing rise to an interaction frequency r'int = 0.35 THz (i.e., 
27rfir'int = 1.4 meV), in excellent agreement with the full 
calculation (underlying color plot). It should be noted 
that this avoided crossing involves a relatively narrow 
mode (the phonon), which can be regarded as a dis¬ 
crete excitation, and a broad one (the bonding plasmon), 
which acts as a continuum, thus leading to Fano inter¬ 
ference [II1II21E1], clearly supported by the presence of 
a deep dip placed in between the two frequencies of the 
resulting hybridized states (see Fig. [^. 


In summary, we here report the observation of plas¬ 
mons in arrays of rings carved in a topological insulator 
( 612803 ) film, as well as their hybridization with a promi¬ 
nent optical phonon of this material in the THz regime. 
As the structures have a thickness that is small compared 
with their lateral dimensions, we are able to describe 


them through classical electrodynamics, using the exper¬ 
imental frequency-dependent surface conductivity of the 
uniform film as input. The latter is taken as the sum of a 
phonon Lorentzian term describing the a 612803 phonon, 
and a Drude term related to Dirac topological surface 
states, with parameters as obtained from a fit to optical 
measurements of unpatterned films. The good agreement 
that we find between theory and experiment through this 
procedure suggests that the response is mainly local. Fur¬ 
thermore, the rings are small compared to the light wave¬ 
length, which allows us to describe them as dipolar scat- 
terers in order to account for the combined effect in the 
arrays. A characteristic bonding/antibonding plasmon- 
mode pattern is observed in the ring response, which is 
strongly affected by the a TI optical phonon, leading to 
Fano profiles in the measured extinction spectra. We 
further observe a strong temperature dependence, we at¬ 
tribute to thermal smearing of the electron distribution 
in the Dirac-fermion states at the insulator surfaces. This 
behavior strongly supports exploiting topological insula¬ 
tor plasmonic resonances in the design of thermally tun¬ 
able terahertz devices. 


METHODS 

Sample fabrication. High quality 6 i 28 e 3 thin-films 
of 60 quintiple layer (QL) thickness (1 QL~1 nm) were 
prepared by a custom-designed 8 VTA M08-V-2 Molec¬ 
ular 6 eam Epitaxy (M 6 E) system using the standard 
two-step growth method [JHlIin]- The lOx lOmm^ AI 2 O 3 
substrates were first cleaned by heating up to 750 °C in 
an oxygen environment to remove organic surface con¬ 
tamination. The substrates were then cooled and kept 
at 110 °C for an initial deposition of 3 quintuple layers 
of 612803 . This was followed by heating to 220°C, at 
which the remainder of the film was deposited to reach 
the targeted thickness. The 8 e and 6 i flux ratio was kept 
constant to approximately 8 e: 6 i= 10 :l in order to mini¬ 
mize 8 e vacancies. Once the films were cooled, they were 
removed from the vacuum chamber and vacuum-sealed 
for transport and handling. 

612803 rings were patterned by electron-beam lithogra¬ 
phy (E 6 L) and subsequent Reactive Ion Etching (RIE). 
The 612803 film was spin-coated with a double layer of 
electron-sensitive resist polymer PMMA (Poly-(methyl 
methacrylate)) up to a total thickness of 1.4 microns. 
The hexagonal ring patterns with different ring sizes were 
then written in the resist by E 6 L. In order to obtain 
a lithographic pattern with re-alignment precision be¬ 
low 10 nm over a sample area suitable for THz spec¬ 
troscopy (^ 10 X 10mm^), we used an electron beam 
writer equipped with a XY interferometric stage (Vistec 
E 6 PG 5000). The patterned resist served as a mask for 
the removal of 612803 by RIE at low microwave power 
of 45 W to prevent heating of the resist mask. 8 ulfur 
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FIG. 2: Plasmon-phonon hybridization in Bi 2 Se 3 TIRs. We represent the extinction maxima (solid symbols) and minima 
(open symbols) of the TIR arrays of[^in a dispersion diagram as a function of inverse outer diameter l/D and light frequency 
V, superimposed on a color plot with the theoretically calculated extinction. The a phonon frequency of Bi 2 Se 3 is shown as a 
horizontal black-dashed line. The two red-dashed curves represent the non-interacting plasmon frequencies calculated without 
including the a phonon (see text and Methods), while the solid curves are obtained from an effective hybridization model (see 
main text). 


hexafluoride (SFg) was used as the active reagent. The 
Bi 2 Se 3 film was etched at a rate of 20 nm/min, and subse¬ 
quently verified by atomic force microscopy (AFM) after 
soaking the sample in acetone to remove PMMA. The 
in-plane edge quality after the RIE process, as inspected 
by AFM, closely followed that of the resist polymer mask 
(i.e., edge roughness smaller than 20 nm). 

Optical characterization. The THz absorption 
spectra were acquired by using a Bruker IFS-125 Michel- 
son interferometer and a liquid-He cooled Infrared-Lab 
bolometer combined with synchrotron radiation at the 
AILES beamline at Synchrotron SOLEIL. The Bi 2 Se 3 
film and a co-planar AI 2 O 3 bare substrate were mounted 
on the cold finger of a pulsed-tube cryoogenerator (Cry- 
omec) and kept at a pressure of about 10“®mbar. The 
extinction coefficient reported in corresponds to 1 — T, 
where T is the film transmittance, defined as the ratio 
between the intensity transmitted through the thin film 
and that transmitted through the bare substrate. 

Modeling the optical response of a topological 
insulator thin film. We intend to simulate the optical 
response of a periodically patterned topological insulator 
film of thickness h small compared with both the lattice 
period a and the lateral dimensions of each structure. 
Under these conditions, the response of the film can be 
described through its effective surface conductivity cr(a;), 
which we assume to be a function of radiation frequency 
oj only. We thus neglect nonlocal effects, which should 


only play a minor role for objects of large lateral dimen¬ 
sions compared with the Fermi wavelength of the electron 
gas in the film. We are describing a spectral region dom¬ 
inated by the presence of an optical phonon, so that its 
bulk dielectric function can be approximated as m 


eti(w) 


^00 




C2 

uj(uj + i 7 ph) ’ 


( 1 ) 


where Wph is the phonon frequency, 7 ph is the phonon 
width, and too is the high-w dielectric function. can 
be equivalently expressed in terms of a bulk conductivity 
as eti = 1 + 47 ricr^°/w, and consequently, it con¬ 
tributes to the 2D surface conductivity of the film a with 
a term ha^^, proportional to the film thickness h. The 
two conducting regions on either side of the film produce 
an additional contribution that should take the form of a 
Drude term at low frequencies. In summary, we can write 
the 2D conductivity of the topological insulator thin film 
as 

, ^ iwe , iwh -S'ph , iwh^_, ^ 

S ca + i7e+4^ u;(w + i7ph)-a;2;^+4^ ^ 

( 2 ) 

where 7 e is the intrinsic decay rate of its surface electronic 
excitations, whereas We is a characteristic frequency that 
should scale with the 1/4 power of the electron gas real 
density, given the Dirac-fermion character of the latter. 
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We note that We should rather be independent of dielec¬ 
tric environment and film thickness h, provided there is 
no overlap between the two metallic regions of the topo¬ 
logical insulator. 

We extract the parameters of the above model for the 
conductivity by means of a fit of the extinction spectra 
of the unpatterned samples, using the standard formula 
for thin films on a substrate m 


[1 + Zoa,{u;)l{n + l)]^ + [Zoa2{io)/{n + 1)]2 ’ 

(3) 

where Zq = 377 is the free-space impedance, cti (w) 
(cr 2 (w)) is the real (imaginary) part of the conductivity, 
and n is the refractive index of the substrate (n = 3.2 for 
AI 2 O 3 in the THz range). The resulting parameters are 
collected in |Tj The effect of these parameters is described 
by the self-consistent optical response elaborated in next 
paragraph. 

Calculation of the transmittance of a periodic 
thin-ring array. We consider thin scatterers of small 
lateral size D arranged in a periodic lattice of period a, 
placed at the interface between materials of refractive 
indices ni = and 712 = For simplicity, we 

focus on normally incident light coming from medium 1 . 
We follow a previously reported formalism to obtain the 
transmission coefficient across this decorated interface as 

where to = 2ni/(ni -|- 712 ) is the transmission coefficient 
for the bare interface, G = 2g/[a^{ei + £ 2 )] + iS” accounts 
for inter-site interactions in the array, g = 5.52 for a 
hexagonal lattice, S = /cA)/{ni+n 2 ), A is the unit¬ 

cell area (i.e., A = (•\/3/2)a^ for a hexagonal lattice), w 
is the light frequency, and a is the polarizability of each 


scatterer including the effect of the interface (i.e., the 
actual induced dipole is given by a times the sum of 
the external field plus the fields scattered by the planar 
interface and the rest of the scatterers). Following an 
electrostatic scaling law that is applicable in the Duj /c 
1 limit, we find [55] 

“ ^ ^ (2T,)/(ei + 62 ) - iujD/a' 

where Aj and Lj are material- and size-independent con¬ 
stants that run over resonant modes j, whereas a is the 
2D conductivity of the thin film (seeabove). This for¬ 
malism is valid under the condition aw/c <C 1 , which is 
fulfilled in this work. 

The extinction is calculated here as 1 — |t/toP where 
to is the transmission coefficient of the sapphire sub¬ 
strate. For simplicity, is given for normal incidence, 
although the calculations here reported are mildly cor¬ 
rected by averaging over a finite range of incidence an¬ 
gles (|0| < sin“^NA for the finite numerical aperture 
NA = 0.11 used in our experiments), using an extension 
of the present theory to oblique incidence PI. 

Parameters for Bi 2 Se 3 rings. We consider rings of 
inner (outer) diameter d (D). These structures are char¬ 
acterized by a bonding mode as well as an antibonding 
mode at higher frequency. In particular, we present cal¬ 
culated results in Fig. |1|2| for hexagonal arrays of rings 
with a D = 3d aspect ratio, characterized by the polar¬ 
izability of with only two terms: [33] Ai = 0.42 and 
Li = 10.3 for the bonding mode; and A 2 = 0.17 and 
L 2 = 47 for the antibonding one. The film thickness is 
h = 60 nm, the lattice period is a = 4d, and the sur¬ 
rounding refractive indices are ni = 1 (air) and 772 = 3.2 
(AI 2 O 3 ). The conductivity of the rings material (Bi 2 Se 3 ) 
is modeled through with temperature-dependent pa¬ 
rameters as shown in Tab. |l|and taking Cao = 30. 
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|a;e/27r (THz) 

7e/27r (THz) 5ph/27r (THz) 

a;ph/27r (THz) 

7ph/27r (THz) 

lOK 

34.1 

1.4 

20.6 

1.9 

0.1 

100 K 

42.3 

2.7 

21.2 

1.9 

0.1 

300 K 

51.0 

4.3 

20.6 

2.0 

0.2 


TABLE I: Fitting parameters extracted to model the conductivity of unpatterned films according to 
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